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ABSTRACT 

As several authors in the past, we calculate optically thick but geometrically thin (and slim) accretion disk models and perform a 
ray-tracing of photons (in the Kerr geometry) to calculate the observed disk spectra. Previously, it was a common practice to ray-trace 
photons assuming that they are emitted from the Kerr geometry equatorial plane, z = 0. We show that the spectra calculated with this 
assumption differ from these calculated under the assumption that photons are emitted from the actual surface of the disc, z = H(r). 
This implies that a knowledge of the location of the thin disks effective photosphere is relevant for calculating the spectra. In this 
paper we investigate, in terms of a simple toy model, a possible influence of the (unknown, and therefore ad hoc assumed) vertical 
dissipation profiles on the vertical structure of the disk and thus on the location of the effective photosphere, and on the observed 
spectra. For disks with moderate and high mass accretion rates (m > 0.01m c ) we find that the photosphere location in the inner, 
radiation pressure dominated, disk region (where most of the radiation comes from) does not depend on the dissipation profile and 
therefore emerging disk spectra are insensitive to the choice of the dissipation function. For lower accretion rates the photosphere 
location depends on the assumed vertical dissipation profile down to the disk inner edge, but the dependence is very weak and thus 
of minor importance. We conclude that the spectra of optically thick accretion disks around black holes should be calculated with 
the ray-tracing from the effective photosphere and that, fortunately, the choice of a particular vertical dissipation profile does not 
substantially influence the calculated spectrum. 
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■ 1. Introduction imation. This simplification leads to a problem: the position 
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of the disk photosphere is calculated only very roughly or is 
There is a consensus that the observed spectra of black hole bina- not calculated at all. In calculating the spectra, one usually 
ries and active galactic nuclei should be explained by accretion assumes that all the emission takes place at the equa torial 
of rotating matter onto black holes. However, the available the- plane . However, in principle ray-tracing routines (e.g. [Bursa 



^ . oretical models of accretion disks do not provide a sufficiently |2006l) should account for the precise location of the emission, 

detailed and accurate description of all physical processes that especially for moderate and high accretion rates. 

are relevant. There is neither a quantitative description of the . . , „ .. „ f , . ,■ 

M F The spectrum computed vertically for each ring (i.e. at 

turbulent dissipation that generates entropy and transports angu- ■ ■• . , x , . , . . . , j. , 

F _ ,f VJ each radius separately), has to be integrated over the disk sur- 

lar momentum, nor a fully self-consistent method to deal with e ^ u ■ ,, c A e , • , u , « „, 

,. . , face. This calls for a need of knowing the global radial 

the radiative transfer in the accreted matter. . e ■• , . . . , u „ 

structure of the accretion disk, consistent with the rings ver- 

Only partial solutions exist. Magnetohydrodynamic tical" structures adopted in the radiative transfer calculations. 

simulations prov ide useful information on the "radial" Expanding along these lines, we are working on constructing 

(lHawlev & Kro lik 2002) and "vertical" jTurneij|2004l) turbulent a new vertical-plus-radial code that inherits the sophisticated 

energy dissipation, but cannot yet simultaneously deal with the treatment of the radiative transfer from the works mentioned 

radiative transfer. On the other hand, all existing methods of above, and at the same time fully incorporates our new radial, 

solving the radiative transfer adopt some simplifying assump- fully relativistic ( in the Kerr geo metry), transonic slim dislQ 

tions. Usually, they divide the flow into separate rings of gas code (described in lSadowskHl2009l) . 

and assume that each ring is in a hydrostatic equilibri um. They 

also assume ad hoc energy dissipation law (e .g. iDavis et aU 

H00l| R6zanska & Madej | 2008t | Idan et alj jjQj). In most cases, > slim disk is a transsonic solution of an optically thick accretion disk 

radiative transfer with all important absorption and scattering described by vertically integrated equations with simplified (limited to 

processes treatment is implemented on the disc surface, while vertical force balance at the surface) treatment of its vertical structure 

the deepest parts of the disk are treated in the diffusion approx- (Abramowicz et al. 1983). 
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Fig. 1 Disk thickness obtained using slim disk model for non- 
rotating BH. Solutions for three mass accretion rates are pre- 
sented: m = 0.001 (dotted line), 0.3 (dashed line) and 0.9/«c 
(solid line). The critical mass accretion rate % corresponds to 
the Eddington luminosity of an accretion disk. 



This paper addresses two questions concerning the effective 
photosphere of optically thick and geometrically thin or slim 
black hole (BH) accretion disk: 

Is a knowledge of the exact photosphere location relevant 
for calculating the disk spectra ? 

Is the photosphere location sensitive to details of dissipative 
processes? 



2. The calculated spectra depend on the location of 
the effective photosphere 

We start with a simple demonstration that the calculated thin 
disk spectra depend on the location of the effective photosphere. 
For this purpose, we use the slim disk soluti ons that have bee n 
recently recalculated in the Kerr geometry bv lSadowskil ([2009). 
We choose three particular models with the accretion rates, 



0.001m c , 0.3m c , 0.9m c . 

Here the critical accretion rate the, defined as 

m c = = 2.23 x 10 18 g ■ s \ 

CK es M Q 



(1) 



(2) 



corresponds to the Eddington luminosity of an accretion disk for 
a non-rotating BH. 

Disk shapes, i.e. the functions z = H(r) describing the ver- 
tical half-thickness, are presented in Fig.Q] For the lowest mass 
accretion rate (O.OOlrac) the H/r ratio for large radii is about 
0.01 while for m = 0.9otc it reaches 0.13. Thus, for the accre- 
tion rates considered here, disks are always geometrically very 
thin, 



H(r) « r. 



(3) 



(Note that for larger accretion rates, the slim disk thick ness could 
be co nsiderably higher.) The models calculated by Sad owskil 
(2009) provide the local flux of radiation F(r,8 e ) in the frame 
of an observer comoving with the disk. In calculating the flux, 
the standard assumption have been adopted: (z) radiation from 
the effective photosphere is locally described by the black-body, 
(if) the flux is limb-darkened by F oul oc 2 + 3 cos 6 e , where 8 e is 
the emission angle relative to the effective photosphere normal 
vector. Using these initial fluxes, we ray-trace photons from the 
effective photosphere to an observer located lOkpc away from 
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Fig. 2 Emerging spectra of accretion disks calculated using the 
slim disk model for three mass accretion rates and different as- 
sumptions about location of the emission for Mbh = 9.4 M . 
The upper panel presents spectral fluxes calculated from the pho- 
tosphere assuming distance to the observer d = lOkpc and the 
inclination angle i = 60°. The ratio of the photon flux obtained 
assuming photosphere emission to the photon flux obtained from 
the equatorial plane is plotted in the lower panel. For the spectral 
fitting purposes energy range 2 -f 20keV is of the major impor- 
tance. 



the central black hole with mass Mbh = 9.4 M (the mass corre- 
sponds to the microquasar 4U 1543-17 while the order of mag- 
nitude of the distance can be considered typical for all BH bina- 
ries) and compute the observed spectra. We take into account all 
relativistic effects in the Kerr geometry (Bursa 2 0061) . 

For the three accretion rates ([1), we calculate the disk spec- 
trum twice: assuming that the effective photosphere coincides 
with the actual disk surface z - H(f), or that it coincides with 
the equatorial plane z — 0, as if the disk would be infinitesimally 
thin. This means that in the first case we start ray-tracing from 
z = H(r), and in the second case from z = 0. 

The resulting spectra are plotted in Fig. [2] There are obvious 
differences between the spectra calculated with the two different 
assumptions about the location of the effective photosphere, i.e. 
the photosphere either at z = H(r), or at z = 0. Even for the 
moderate mass accretion rate 0.3mc, the "z = 0" spectral flux 
differs from the "z = H{r)" flux by about 20% at 20keV, and 
for higher accretion rates the difference is much higher. Such 
differences should certainly be considered as highly significant. 

Thus, our result demonstrates that knowing the precise loca- 
tion of the effective photosphere is necessary in calculating the 
BH slim disk spectra. The effective photosphere lies certainly 
somewhere between the z = H(r) and z — surfaces. In princi- 
ple, its location could strongly depend on details of the vertical 
dissipation, and this would be a rather bad news — these de- 
tails are still not sufficiently well known and therefore ad hoc 
assumed in models. 



A. Sadowski et al.: Accretion disk photosphere 



3 



In the rest of this paper we argue that such a situation is un- 
likely. Using a simple toy model for the vertical structure of ge- 
ometrically thin, optically thick accretion disk, we show that the 
photosphere location is not highly sensitive to dissipation. Note, 
that for very thin © stationary disks the same is true for the total 
flux F(r) emitted from a particular radial location — it doe s not 
depend on dissipation processes (Shakura & Sunvaev 1973b . 

These results strengthen one's confidence in estimating the 
black hole spin by fitting the observed spectra o f black hole 
sourc e s to these calculated theoretically (see e.g. IShafee et"aT1 
2006; iMiddleton etail 120061; iGou et al.1 l2009t) . However, they 
also indicate that the fitting procedures should be improved to 
include ray-tracing from the actual location of the photosphere. 



3. A simple model of the vertical structure 

Our simple model assumes that the accretion disk is geometri- 
cally thin (01 and that one may consider radial and vertical disk 
structure separately. 



3.1. Radial equations 

We do not solve radial equations, assuming instead that the rota- 
tion is strictly Keplerian (in the Kerr geometry), 



O = Q 



K 



r 3/2 . 



( y/GM/c 



GM 1 



(4) 



and that the flux F{r) follows from the mass, energy and angular 
momentum conservation (does not de pend on radial dissipation) 
and is given by the s tandard formula dNovikov & Thorn"g|l973 ; 
Page & . Thornelll974h . 



_, . 3 GMM 
F(r) = r- 



Q 



SC*' 2 ' 



(5) 



Herein M is the black hole mass, a is its specific angular momen- 
tum, and ${(r, a), S(r, a), C(r, a), V(r, a), 8(r, a), Q{r, a) are rel- 
ativistic corr ection factors, defined as explicit functions of their 
arguments in lPage & Thornd d!974h . 

3.2. Vertical equations 

We describe the vertical structure of the BH accretion disk in 
optically thick regime by the following equations: 
(i) The hydrostatic equilibrium: 

dP , 

— = -p£l\@z (6) 
dz 

where P is the sum of the gas and radiation pressures: 

P = Pgas + Prad = kpT + jClT 4 (7) 

and £l 2 K @z being the vertical component of the gravitational 
force of the central object calculated using the relativistic cor- 
rection factor: 



~~ W 2 C ' 

(ii) The energy generation equation: 



dF 
dF 



tits 



(8) 



(9) 



where F is the flux of energy generated inside the disk at a rate 
given by Qm s - the dissipation rate which will be discussed in 



(iii) The generated flux of energy is transported outward 
through diffusion of radiation and/or convection according to the 
following thermodynamic al relation: 



dlnr 
dlnP 



(10) 



where V is the thermodynamical gradient which can be either 
radiative or convective: 



V = 



V rad dla V rad < V ad 



dla V rad > V, 



ad 



(ii) 



The radiative gradient V rad is calculated using diffusive approx- 
imation: 



^ 1 rod = 



3k r PF 
16o-T 4 Cl 2 K gz 



(12) 



where k r is the mean Rosseland opacity (see 83.31 ) and <x is the 
Stefan-Boltzmann radiation constant. 

(iii) When the temperature gradient is steep enough to ex- 
ceed the value of the adiabatic gradient V ad we have to consider 
the convective energy flux. The convective gradient V conv is cal- 
culate d using the mixing length theory introduced by Paczvhski 
d!969t) . Herein we take the following mixing length: 

H ml = \.M P (13) 

with pressure scale height Hp defined as (lHameurvet all 1998): 

Hp = — ^= — (14) 



pQ. 2 K z + yfPpQ. K 
The convective gradient is defined by the following formula: 
^ ' conv — * conv + (V rarf - V ad )y(y + V) (15) 
where y is the solution of the equation: 
9 



y 3 + y y 2 + y2 y . 



V = 



(16) 



4 3 + t 2 , 

ml 

with the typical optical depth for convection r m / = Kp.pH m \ and 
V given by: 

Wrad ~ V fld ) (17) 



1 



( 3 +t 2 ^ 



3t„ 



512o- 2 T 6 H P \31nT 



The thermodynamical quantities Cp, V ad and (d In p/d In T)p are 
calculated using standard formulae (e.g. IChandrase kharl Il967l) 
assuming solar abundances (X = 0.70, Y = 0.28) and taking into 
account, when necessary, the effect of partial ionization of gas 
on the gas mean molecular weight. 

(iv) To close the set of equations describing the vertical struc- 
ture of an accretion disk we have to provide boundary condi- 
tions. At the equatorial plane (z = 0) we put F = while at the 
disk photosphere we r equire F = crT 4 with the flux F given by 
lPage&Thomel(ll974l) . 

3.3. Opacities 

In this work we consider optically thick accretion disks. 
Therefore, we use Rosseland mean opacities kr. The opaci- 
ties as a function of d ensit y and tempe r ature are taken from 
lAlexander et al.l d 1983b and ISeaton et alj d!994l) . For tempera- 
tures and densities out of both domains we interpolate opacities 
between these two tables for intermediate values. 
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4. Numerical method 

The set of ordinary differential equations defined in $3.21 is 
solved using Runge-Kutta method of the 4th order. We start inte- 
grating from the equatorial plane (z - 0) where we assume some 
arbitrary central temperature Tq, density pc and set F — 0. We 
integrate given set of equations until we reach the photosphere 
which is defined as a layer where the following equation (corre- 
sponding to the optical depth t = 2/3): 



k r P = -£l K @z 



(18) 



is satisfied. If the temperature T and flux F does not satisfy the 
outer boundary condition (F = oT 4 ) we tune up the value of 
the assumed central density pc and again integrate starting from 
the equatorial plane. The converged solution is usually obtained 
after a few iteration steps. To get the photosphere profile h(r) 
we look for the disk solutions at a number of radii in the range 
r ms < r < 200M with r ms being the radius of the marginally 
stable orbit. At each radius we additionally search for the value 
of the central temperature which determines the solution with 
the emitted flux equal to the Novikov & Thorne value. 



5. The photosphere location for a few ad hoc 
assumed vertical dissipation profiles 

Herein we apply our simplified model of vertical structure (|3) 
to calculate photosphere location for different dissipation pre- 
scriptions. We test four families of dissipation functions de- 
scribed in the following section. 



5.1. Dissipation profiles 

First of all we apply the standard aP prescription given by 
Shak ura & Sunvae"vl(ll973l) . The authors assumed that the source 
of viscosity in accretion disk is connected with turbulence in gas- 
dynamical flow the nature of which was unknown at that time. 
They based on the following expression for the kinematic vis- 
cosity coefficient v: 



V ~ Vturbhurb 



(19) 



(where v tU rb an d hurb stand for typical turbulent motion velocity 
and length scale, respectively). Shakura & Sunyaev assumed that 
the velocity of turbulent elements cannot exceed the speed of 
sound c s as well as their typical size cannot be larger that the 
disk thickness h. Taking into account the expression for the t r $ 
component of the viscous stress tensor in a Newtonian accretion 
disk: 



fr0 = --jPVtl 

and the standard form of the vertical equilibrium: 



,1/2 



r 

GM 



1/2 



one can limit the t r<t) in the following way: 
- t r<t> < pQ.c s h x pc 2 s K P 



(20) 



(21) 



(22) 



Therefore, we may introduce a dimensionless viscosity parame- 
ter a satisfying the condition a < 1 : 



This expression for viscosity is called the a prescription and has 
been widely and successfully used for many years in accretion 
disk theory. In this case the flux generation rate dF/dz (Eq. [9j, 
given in the local rest frame by: 



dF _ 3DS 

-j- - 0> )W f (r)(0) - -^-^-t (rm Ll K 

is expressed as: 

dF 3 DS 
Qdis = -7- = - —aP£l K 
dz 2 C 



(24) 



(25) 



work we apply the regular 
a few arbitrary values of the 



aP prescrip- 
a parameter: 



In this 
tion for 
a = 0.001,0.01,0.1,0.2,0.5, 1.0. 

In the last decade several authors (e.g. ISincell & Krolikl 

[19981: iHubenv et al.l l200Tt iDavis et all 120051 iHui et alJ 12005) 
used to apply a constant dissipation rate per unit mass. We fol- 
low them investigating a class of models using the following for- 
mula for the flux generation rate: 



Qdis - 



dF _ F 
dz ni(/ 



(26) 



Due to the fact that this expression requires not only the total 
emitted flux Fq but also the total surface density mo one has to 
provide the latter basing on other disk solutions. In our work we 
use models calculated using regular a prescription for different 
values of a. 

Recent numerical MHD simulatio ns of stratified acc retion 
disks in shearing box approximation dDavis et al. 2009) show 
that dissipation rate can be well approximated by the following 
formula: 



Qdis - 



dF 1 F {) 1 



dz 2 s t /m ~ sjm 



(27) 



We apply it in another class of models (mT l l 2 profile). Similarly 
like in the constant dissipation rate case we impose surface den- 
sity profile obtained using aP models. 

Several authors investigated accretion disks with heat gener- 
ation concen trated near the equatorial plane or close to the dis k 
surface (e.g. lPaczvnski & Abramowicz 1982b IPaczvn ski 1980). 
Therefore, we implement one more class of dissipation prescrip- 
tions with arbitrary positioned dissipation maximum. We use 
the following expression 



Qdis - 



dF 
dF 



dF 
dF 



3D® 

= aP£l K e 

2 C 



(28) 



being flux generation rate for aP models re-scaled with Gaussian 
function centered at zo with dispersion parameter d given by: 



zo/r = d tl 



d/r = d di! 



sp 



(29) 



t r <b = -aP 



(23) 



where d max and ddisp are constant for a given disk model. 

In Fig. [3] we present dissipation profiles for a few exem- 
plary models. The upper panel presents flux generation rate ver- 
sus the vertical coordinate. The aP model (solid line) has bell- 
shaped dissipation profile where most of energy is generated at 
z < 0.5//. The model with constant dissipation rate per unit mass 
(thin dashed line) exhibits very similar behavior. The dissipation 
in no longer concentrated around the equatorial plane for the 
m ^ 2 model (thick dashed line) - the maximal flux generation 
is located at z — 0.4//. However, the dissipation profile is still 
extended through the whole disk thickness. The arbitrary posi- 
tioned models (dot-dashed lines) show opposite behavior: the 
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Fig. 3 Dissipation profiles for a few representative models at 
r - 2QM. The upper panel presents dissipation functions versus 
the vertical coordinate z while the bottom panel versus surface 
density m(z) - J H p(z)dz'. On both panels left and right edges 
correspond to disk surface and equatorial plane, respectively. 
The models represent four families: aP (A3, solid line), constant 
dissipation rate per unit mass (B2, thin dashed line), mT 112 pro- 
file (C2, thick dashed line) and arbitrary positioned dissipation 
function (Dl and D5, dash-dotted lines). 



energy generation is concentrated around an arbitrary given lo- 
cation. In case of the models presented in Fig.[3]the dissipation 
is confined to altitudes close to the equatorial plane (z = 0) and 
z = 0.5//. 

Another point of view on the dissipation profiles is presented 
in the bottom panel of Fig. [3] The profiles are plotted against the 
fractional mass surface density m/mo. The right edge of the fig- 
ure (m = mo) corresponds to the equatorial plane. Two families 
of models, with constant and proportional to m~'/ 2 dissipation 
profiles, are represented by straight lines (they are power-law 
functions of m). As on the upper plot, one can notice similar- 
ity between the aP and constant dissipation rate models. The 
maximum of the flux generation rate for the arbitrary positioned 
model centered at z — 0.5// corresponds to the mass depth 
m = O.lmo. Obviously, the model centered at the equatorial 
plane has maximal dissipation rate at m = mo- 



Details of model parameters are given in Tab. [T] The location of 
the disk photospheres is presented in Fig. [4] The first four pan- 
els present results for four classes of dissipation profiles in case 
of a non rotating BH. Common behaviour is clearly visible: all 
photosphere locations coincide at radii lower than 30M. Outside 
this radius the location of the photosphere depends on the as- 
sumed dissipation profile and vary between H/r as 0.03 for the 
aP model with the highest value of a — 1 .0 and H/r as 0.06 for 
the lowest value a = 0.001. It is of major importance to under- 
stand why all photosphere profiles coincide in the inner region 
which corresponds to the radiation pressure and electron scatter- 
ing dominated regime of an accretion disk. The location of the 
photosphere is defined by Eq. [18] Wherever the radiation pres- 
sure and electron scattering dominate in a disk the left hand side 
of that formula depends only on the temperature which is con- 
nected at the photosphere with the outcoming flux by F = crT 4 . 
Therefore, the formula for the energy flux (Eq. [5) determines 
uniquely the location of the photosphere in the inner region of 
a disk independently of the assumed dissipation function. Under 
these assumptions we get from Eq. QjU 

H = ^ (30) 

which describes perfectly the photosphere profiles presented in 
Fig.0forr < 30M. 

Disk thickness profiles for a case of a rotating BH for three 
aP models are presented in the bottommost panel. They exhibit 
very similar behaviour: the photosphere location does not de- 
pend on the viscosity prescription for small radii corresponding 
to radiation pressure and electron scattering dominated region 
which, in case of a rotating BH with a* = 0.9, extends inside 
r = 40M. Therefore, we may infer that such behaviour is general 
in accretion disks which exhibit radiation pressure dominated 
regimes. This statement is of paramount importance due to the 
fact that most of the emission from an accretion disk take place 
(in non-rotating case) between 6 and 2QM. As we have proven, 
for such radii the photosphere location should not depend on as- 
sumed dissipation profile and, therefore, the spectrum is likely to 
be independent of dissipation profile, either. However, one can- 
not make sure without performing full ray-tracing which would 
account for flux emitted from dissipation dependent regions cor- 
rectly. Results o f this procedure will be disc ussed in the follow- 
ing section. As Shakura & Sunvaev ( 1973) have shown for the 
lowest mass accretion rates the gas dominated region of an ac- 
cretion disk may extend all the way down to the inner edge of a 
disk. For such case the behavior described in the previous para- 
graph is not expected: the photosphere profile should depend on 
the dissipation function at all radii. To account for this fact we 
compare photosphere profiles of accretion disks with low mass 
accretion rates (Fig. [5j. For m - Q.lthc the radiation dominated 
region extends upto r ~ 40M. For lower accretion rate (0.05mc) 
it is confined to r < 20M, while for the lowest (O.Olrac) it never 
appears. In that case the photosphere location indeed depends 
on the assumed viscosity prescription even for radii lower than 
2QM where most of the emission comes out. However, one has 
to bear in mind that H/ r ratio for such low mass accretion rate 
is even below 0.015 and the impact of the photosphere location 
on the emerging spectrum should be insubstantial. It is studied 
in details in the following section. 



5.2. The photosphere location 

We test a number of models representing all four classes of dis- 
sipation profiles for a moderate mass accretion rate m = O.lmc. 



5.3. Spectra 

To assess the importance of different vertical dissipation pro- 
files on emerging spectra we calculate them basing on photo- 
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Table 1 Model Assumptions form = 0.1 the 



Model 


Model family" 


a 


dmax 




a' 


Al 


aP 


0.001 


- 


- 


0.0 


A2 


aP 


0.01 






0.0 


A3 


aP 


0.1 






0.0 


A4 


aP 


0.2 






0.0 


Aj 


aP 


U.J 






A A 

U.U 


Ao 


aP 


1 A 
1 ,U 






A A 

U.U 


AO 1 


tv P 
{it 


U.Ul 






yl.J 


A*\9 


rv P 


1 






Q 

V/. 7 


AS 3 


aP 


0.5 






0.9 


Bl 


pnnctQnt Hiec raff 1 


01 b 

U.U 1 






o n 


R9 

D _ 


t-UllMclllL Ulaa. IdLC 


U. 1 









DT 
DJ 


constant diss, rate 


U.J 






A A 

u.u 


L 1 


*™-t/2 


U.Ul 






A A 

U.U 




in 


n 1 b 

U. 1 






A A 
U.U 


C3 


-1/2 


0.5 b 






0.0 


Dl 


arbitrary positioned 


0.1 


0.0 


0.001 


0.0 


D2 


arbitrary positioned 


0.1 


0.02 


0.001 


0.0 


D3 


arbitrary positioned 


1.0 


0.0 


0.001 


0.0 


D4 


arbitrary positioned 


1.0 


0.02 


0.001 


0.0 


D5 


arbitrary positioned 


1.0 


0.02 


0.0005 


0.0 


D6 


arbitrary positioned 


1.0 


0.02 


0.005 


0.0 



"Details of model assumptions are given in 35.11 

'Calculated basing on surface density profiles obtained in the aP model 
with given a and a' 

sphere profiles obtained using our model. For comparison, we 
also account for spectra calculated assuming emission from the 
equatorial plane. Spectra for models with mass accretion rate 
th - O.lthc are presented in Fig. [6] The upper panel presents 
total spectra for non- and rapidly-rotating BH. The spectral pro- 
file for a* = 0.9 case extends to higher energies due to the fact 
that accretion disk inner edge moves closer to the central ob- 
ject with increasing BH angular momentum extracting more en- 
ergetic photons. Spectra for all the photosphere profiles almost 
coincide and differ significantly from spectra calculated assum- 
ing emission from the equatorial plane. It is clearly visible on 
the lower panel where we plot ratio of spectral fluxes emitted 
from the photosphere to emitted from the equatorial plane for 
two most extreme photosphere profiles (compare Fig. |4j for each 
value of BH angular momentum. The spectral profiles are iden- 
tical up to a factor of 1%. The agreement is the best for ener- 
gies of few keV. One may conclude that for accretion disks with 
sufficiently high accretion rate to form the inner radiation pres- 
sure supported region, the spectra depend only insignificantly 
on vertical dissipation p rofile in the mo st interesting range of 
energies (2.5 - 20keV; iGou et alJl2009h . As was discussed in 
35.21 for the lowest mass accretion rates (< 0.05toc) the inner, 
radiation pressure dominated disk region never appears and the 
photosphere location depends on the dissipation profile down to 
the disk inner edge. As it has been pointed out, such regime of 
accretion rates implies small values of the H/R ratio. To assess 
the impact of the dissipation profiles on the spectra we plotted 
ratios of corresponding spectral profiles for low and very low 
mass accretion rates (Fig. [7]). The upper panel presents the ra- 
tios as observed by an observer perpendicular to the disk plane. 
The ratios approach 1 with decreasing mass accretion rate. The 
convergence is very slow due to the fact that disk thickness in 
disk outer region (gas pressure and free-free scaterings domi- 
nated) depends weakly on the mass accretion rate (H ~ m 3 ^ 20 , 
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Fig. 4 Photosphere profiles for different model families for th = 
O.lthc- The uppermost panel presents profiles for the aP models. 
The second one for models with constant dissipation rate per unit 
mass. The third one for m ^ 2 dissipation profile. The forth is for 
arbitrary positioned dissipation while the bottommost is for aP 
models in case of a rotating BH {a* = 0.9). BH mass was 9.4 M©. 
Model parameters are given in Table Q] 

Shak ura & Sunvae"vlll973l) . The difference in spectral fluxes be- 
tween models with th = O.lmc and O.OOOlmc is no larger than 
0.1% at 20keV. For inclined observer (the bottom panel) the de- 
partures from the equatorial plane model are much more signif- 
icant (up to 10% at 20keV even for the lowest mass accretion 
rate). However, the difference between the spectral fluxes for the 
models with the two lowest mass accretion rates (corresponding 
to H/R = 0.007 and 0.015, respectively) is still of the order of 
1% or even smaller. We conclude that in the lowest accretion 
rate regime changes of the photosphere location caused either 
by different mass accretion rates or different dissipation profiles 
are insignificant. However, one has to keep in mind that even for 
mass accretion rates as small as O.OOOlmc taking into account 
the non-zero photosphere location is important for inclined ob- 
servers as the spectral flux at high energies approaches the equa- 
torial plane limit very slowly. 

6. Discussion 

Results presented here are very encouraging in the context of fit- 
ting the calculated optically thick spectra of slim accretion BH 
disk to the observed spectra. Our results imply that the major un- 
certainties of the accretion disk models, in particular the vertical 
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Fig. 5 Photosphere profiles for aP models for three different 
mass accretion rates (m - 0.10, 0.05 and 0.01;«c) and two val- 
ues of the a parameter (0.1 and 1 .0). The regions where the pho- 
tosphere profiles for a given accretion rate coincide are radiation 
pressure dominated. BH was non-rotating with M = 9.4 M e , 
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Fig. 6 Emerging spectra calculated for different models. The up- 
per panel presents number of photons per energy crossing every 
second a unit surface assuming d - lOkpc. Spectra calculated 
assuming emission from the equatorial plane (z em = 0) as well 
as from the photosphere surface (z em = H) are presented for 
non-rotating (a* = 0) and highly-spinning (a* - 0.9) BH. The 
spectra for different dissipation profiles are indistinguishable - 
they all coincide with dot-dashed and dotted lines (depending on 
BH spin). The bottom panel presents ratio of the spectral flux 
(defined as above) calculated from the photosphere to the flux 
calculated from the equatorial plane for the models with extreme 
photosphere profiles (see j j5.2| i. The inclination angle i = 60°. 



dissipation profiles, have a rather small influence on the calcu- 
lated spectra of steady disks with low accretion rates. This is a 
good news for these who use spectral fitting to estimate the black 
hole spin. 

One should have in mind, however, that there are several ef- 
fects that have been not taken into account in our simple model, 
and more theoretical work is needed before fully believable 
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Fig. 7 Ratios of spectral profiles (photosphere to equatorial) for 
models with low mass accretion rates for the perpendicular (i = 
0°, upper panel) and i = 60° (lower panel) observers. 



methods of calculating slim disk spectra will be at hand. One 
of the most obvious (and also quite challenging) theoretical de- 
velopments needed is a self-consistent and simultaneous solving 
the vertical and radial structures of slim, optically thick black 
hole accretion disk. We are working on this problem. 
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